Direct Fisher inference of the quartic oscillator's eigenvalues 
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Abstract 

It is well known that a suggestive connection links Schrodinger's equation (SE) and the information- 
optimizing principle based on Fisher's information measure (FIM). It has been shown that this 
entails the existence of a Legendre transform structure underlying the SE. Such a structure leads 
to a first order partial differential equation (PDE) for the SE's eigenvalues from which a com- 
plete solution for them can be obtained. As an application we deal with the quantum theory of 
anharmonic oscillators, a long-standing problem that has received intense attention motivated by 
problems in quantum field theory and molecular physics. By appeal to the Cramer Rao bound we 
are able to Fisher-infer the particular PDE-solution that yields the eigenvalues without explicitly 
solving Schrodinger's equation. Remarkably enough, and in contrast with standard variational 
approaches, our present procedure does not involve free fitting parameters. 
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1. INTRODUCTION 



It is well-known that a strong link exists between Fisher' information measure (FIM) / [lj and 
Schrodinger wave equation (SE) [2|-l8|. In a nutshell, this connection is based upon the fact 
that a constrained Fisher-minimization leads to a SE-like equation [l-8]. In turn, this implies 
the existence of intriguing relationships between various SE-facets, and Jaynes's maximum 
entropy principle. In particular, basic SE-consequences such as the Hellmann-Feynman and 
the Virial theorems can be re-interpreted in terms of a special kind of reciprocity relations 
between relevant physical quantities, similar to the ones exhibited by the thermodynam- 
ics' formalism via its Legendre-invariance property This fact demonstrates that a 
Legendre-transform structure underlies the non-relativistic Schrodinger equation. As a con- 
sequence, the possible energy-eigenvalues are now seen to be constrained by such structure 
in a rather unsuspected way js — 8 1 , a fact that allows one to obtain a first-order differential 
equation, unrelated to Schroedinger's equation {?], that energy eigenvalues must neces- 
sarily satisfy. The predictive power of this new equation will be explored here. 
We will apply our formalism here to the quantum anharmonic oscillator, which is the paradig- 
matic testing-ground for new approaches to Schroedinger eigenvalue equation. Besides their 
intrinsic conceptual and mathematical interest, anharmonic oscillators have received consid- 
erable attention over the years due to their practical relevance in connection with several 
areas of physics, such as quantum field theory and molecular physics, among others. In 
this kind of systems, the most intense focus has been traditionally concentrated upon the 
quartic oscillator. General accounts containing illuminating references on this problem may 
be found, for instance, in ^, [loj]. Note that a perturbation series solution to this problem in 
powers of the anharmonicity-parameter A is divergent-asymptotic for all A > Specif- 
ically, we will apply our procedure to treat the quartic anharmonic oscillator. By appeal to 
the Cramer Rao bound we obtain the particular solution that leads to the system's eigen- 
values without need of explicitly solving Schrodinger's equation. More importantly, we do 
not need to app eal to any arbitrary or empirical parameter, as is common practice in other 



treatments 



12|. 
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2. BASIC IDEAS 



Let x be a stochastic variable and f(x,9) the probability density function (PDF) for this 
variable, which depends on the parameter 9. If an observer were to make a measurement of 
x and had to best infer 9 from such measurement, calling the resulting estimate 9 = 9(x), 
one might ask how well could 9 be determined. Estimation theory [l| tells us that the best 
possible estimator 9 (x), after a very large number of x-samples is examined, suffers a mean- 
square error Ax from 9 obeying the rule / (Ax) 2 = 1, where / is an information quantifier 
called the Fisher information measure (FIM), a non linear functional of the PDF that reads 

/ = J dxf(x,9)^ ln[/M)]j . (1) 

Any other estimator must have a larger mean-square error (all estimators must be unbiased, 
i.e., satisfy (#(x)) = 9). Thus, FIM has a lower bound. No matter what the parameter 9 
might be, / has to obey 

I {Ax) 2 > 1, (2) 

the celebrated Cramer-Rao bound [lj]. 

In the case of physical Fisher applications the particular instance of translational families 
merits special consideration. These are mono-parametric distribution families of the form 
f(x, 9) = f(x — 9), known up to the shift parameter 9. All family members exhibit identical 
shape. For such families we get 

x)y 

dx 



!(x) ( d ±m ]ix . (3 ) 



Focus attention now a system that is specified by a set of M physical parameters [i k . We 
can write \i k = (A k ), with A k = A k (x). The set of fi k - values is to be regarded as our prior 
knowledge. It represents our available empirical information. Let the pertinent probability 
distribution function (PDF) be f(x). Then, 



(A k ) = J dx A k (x) f(x), k — l,...,M. (4) 

In this context it can be shown (see for example {2, 3|) that the physically relevant PDF 
f(x) minimizes FIM subject to the prior conditions and the normalization condition. Nor- 
malization entails / dxf(x) = 1, and, consequently, our Fisher-based extremization problem 
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adopts the appearance 



5 ^1 — a J dx f(x) ~^2^kJ dx A], 



(x) f(x) 



0. 



(5) 



where we have introduced the (M + 1) Lagrange multipliers A& (A = a). In Ref. [2] on 
can find the details of how to go from fl5]) to a Schrodinger's equation (SE) that yields the 
desired PDF in terms of the amplitude ip(x) defined by f(x) = ip(x) 2 . This SE is of the 
form 



1 d 2 



2 dx' 



+ U(x) 



TP 



a 

77 ^ 



U(x) 



1 



M 



- X k A k( X ), 



(6) 



k=l 



and can be formally interpreted as the (real) Schrodinger equation (SE) for a particle o" 
unit mass (h = 1) moving in the effective, "information-related pseudo-potential" U(x) [2] 
in which the normalization-Lagrange multiplier (a/8) plays the role of an energy eigenvalue. 
The Xk are fixed, of course, by recourse to the available prior information. In the case of 
one-dimensional scenarios, ip(x) is real 13) and 



so that using the SE 



1 dln^p 
y dx J 

we obtain 



M 



a. 



(8) 



fc=i 



Legendre structure 



The connection between the variational solutions / and thermodynamics was established in 
Refs. Q and Q in the guise of reciprocity relations that express the Legendre-transform 
structure of thermodynamics. They constitute its essential formal ingredient Ij| and were 
re-derived a la Fisher in j^J by recasting (jSj) in a fashion that emphasizes the role of the 
relevant independent variables, 



I ((A,),..., (A M )) = « + E X k (A k ). 



(9) 



k=l 



Obviously, the Legendre transform main goal is that of changing the identity of our relevant 
variables. As for I we have 

M 



a(Ai, . . . , A M ) = 1 ~J2 ( A k) 
k=i 



(10) 



Q: 



so that we encounter the three reciprocity relations (proved in [2j) 

da - (A )■ 91 -X ■ 81 -VA 9{Ak) (11) 

the last one being a generalized Fisher-Euler theorem. 



3. FISHER MEASURE AND QUANTUM MECHANICAL CONNECTION 

The potential function U(x) belongs to £2 and thus admits of a series expansion in the basis 
x, x 2 x 3 , etc. 15] . The A k {x) themselves belong to £2 as well and can also be series-expanded 
in similar fashion. This enables us to base our future considerations on the assumption that 
the a priori knowledge refers to moments x k of the independent variable, i.e., (Ak) = (x k ), 
and that one possesses information about M of these moments (x h ). Our "information" 
potential U then reads 



U(x)=-^\ k x k . (12) 

k 

We will assume that the first M terms of the above series yield a satisfactory representation 
of U(x). Consequently, the Lagrange multipliers are identified with U(x)'s series-expansion's 
coefficients. 

In this Schrodinger-scenario the virial theorem states that jsj 
I B 2 \ IB \ 1 M 

s =- '^"riS'^' (13) 



and thus, from (J7|) and 
sure can be arrived at [5] 



131) a useful, virial-related expression for Fisher's information mea- 



i = - £2 Afc(xfc) ' (14) 

k=l 

Also, substituting the above /-expression into (jSJ) and solving for a, we obtain 

M 



a = — 

k 



M ( k\ 

E [l + ^jXk(x k ). (15) 



a (I) is explicit function of the M Lagrange multipliers - £7(x)'s series-expansion coefficients 
Afc (associated to the physical parameters (x k )). Eqs. ( !T4|) and (|T5|) encode the information 
provided by the virial theorem 



Fisher-Schroedinger Legendre structure 



Interestingly enough, the reciprocity relations (RR) ( ITTj) can be re-derived on a strictly pure 
quantum mechanical basis [5], starting from 

1. the quantum Virial theorem [which leads to Eqs. (fl^|) and ( TT5|) ] plus 

2. information provided by the quantum Hellmann-Feynman theorem. 

This fact indicates that a Legendre structure underlays the one-dimensional Schroedinger 
equation [5]. Thus, with (A k ) = (x k ), our "new" reciprocity relations are given by 

£ = w^ = ^ w = < 16 > 

oX k o{x k ) d\i Y dX i 

FIM expresses a relation between the independent variables or control variables (the prior in- 
formation) and I. Such information is encoded into the functional form / = I((x l ), (x M )). 
For later convenience, we will also denote such a relation or encoding as {J, (x k )}. We 
see that the Legendre transform FIM-structure involves eigenvalues of the "information- 
Hamiltonian" and Lagrange multipliers. Information is encoded in / via these Lagrange 
multipliers, i.e., a = a(X±, ...\m), together with a bijection {I, (x k )} < — > {a, Xk}- 

Two scenarios 

In a j/, - scenario, the are functions dependent on the (x fc )-values. As shown 



in jg], substituting the RR given by (jT6l) in pl|) one is led to a linear, partial differential 
equations (PDE) for /, 

Xk = W) ~* I = ~h2 {x) W)- (17) 

and a complete solution is given by 

I((x'),...,(x M ))= f C k \{x k )[ V \ (18) 

k=l 

where C k are positive real numbers (integration constants). The I - domain is Dj = 
{((x 1 ),..., (x M ))/(x k ) e Ko}- Eq. (HHD states that for (x k ) > 0, J is a monotonically 
decreasing function of (x k ), and as one expects from a "good" information measure I is 
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a convex function. We may obtain Xk from the reciprocity relations ( TT6|) . For (x ) > 

one gets, 

Xt = W) = -l c - {xt) ~ i2+t '" k < • (19) 

and then, using (jH}, we obtain the a - normalization Lagrange multiplier. 
The general solution for the I - PDE does exist and its uniqueness has been demonstrated 



via an analysis of the associated Cauchy problem [6]. Thus, Eq. (118p implies what seems to 
be a kind of "universal" prescription, a linear PDE that any variationally (with constraints) 
obtained FIM must necessarily comply with. 

In the {a, Xk} scenario, the (x k ) are functions that depend on the A^-values. As we showed 
in [7], an analog a-PDE exists. Substituting the RR given by (IT5|) in (fT5j) we are led to 

and a complete solution is given by 

M 

a(A 1 ,...,A A/ )= £ D fc |A fc | 2/(2+fc) , (21) 
fc=i 

where the D^s are positive real numbers (integration constants). The a-domain is D a = 
{(Ai, • • • , A^)/Afe G 3?} = 3? M . Also, Eq.( l2~T|) states that for A^ < 0, a is a monotonically 
decreasing function of the Xk, and as one expect from the Legendre transform of /, we end 
up with a concave function. We may obtain the (x fc )'s from the reciprocity relations (I16p . 
For Xk < one gets 

<**> = - k = (^i) c ' i^i- t/,2+t) > ' 22 > 

and then, using ([8]) one us able to build up /. 

The general solution for a - PDE exists. Uniqueness is, again, proved from an analysis of 
the associated Cauchy problem Q . Thus, Eq. (TPS]) implies once more a kind of "universal" 
prescription, a linear PDE that all SE-eigenvalues must necessarily comply with. 

The mathematical structure of the Legendre transform leads to a relation between the in- 
tegration constants Ck and Dk pertaining to the / and a expressions, respectively, given 
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by (ITS]) and (12TI) . In [3] we studied with some detail this relation. In our two scenarios, 
{/, and {a, A^}, we have [7] 



Cu 



k c 



with D k (2+h) = C k k =Fl (23) 



Consequently, expressions flTHj) and ( l2~Tj) take the form, 



M u 

£ I 

fc=i z 



2/fc 



E ^ [*i |A fe |] 2/(2+fc) . 



(24) 



fc=i 



The reciprocity relations ( 1191) and ( 1221) can thus be economically summarized in the fashion 



F k 2 



(25) 



4. PRESENT RESULTS 

The reference quantities 

The essential FIM feature is undoubtedly its being an estimation measure known to obey 
the Cramer Rao (CR) bound of Eq. ()2]) [lj. Accordingly, since our partial differential 
equation has multiple solutions, it is natural to follow Jaynes's MaxEnt ideas and select 
amongst them the one that optimizes the CR bound, that constitutes the informational 
operative constraint in Fisher's instance. Of course, Jaynes needs to maximize the entropy 
instead. We will also, without loss of generality, renormalize the reference quantities F k . This 
procedure is convenient because it allows us to regard these quantities as statistical weights 
that optimize the CR-bound. In other words, our procedure entails that we extremize 

M k 

f(F 1 ,-..,F M ) = I ((x 2 ) - (x) 2 ) = E | 

k=l 

with the constraint 

M 

<f>(F ir ..,F M ) = £ F 2 ' k = l. (27) 

k=l 

We are going to apply now the preceding considerations so as to obtain the eigenvalues of 
the quartic anharmonic oscillator. 



F k 

K**>l 



((x 2 ) - (x) 2 ) 



(26) 



S 



Quartic anharmonic oscillator 



The Schrodinger equation for a particle of unit mass in a quartic anharmonic potential reads, 

ip = E ip. (28) 



id 2 i, 2 1 4 

2 kX + 2 XX 



where A is the anharmonicity constant. According to [7|, |8J|, we can ascribe to ( 12 8 D a Fisher 
measure and make then the following identifications: a = 8E , A2 = —4 k , A4 = —4 A. 
Accordingly, we have, in the {a, Xk} - scenario [Cf. (124"])]. 



a: 



2F 2 1/2 |A 2 | 1/2 + 3F 4 1/3 |A 4 | 1/3 



(29) 



The functions / and (ft defined by (T2"6"]) and (I2T]) . respectively, can here be recast [using (123]) ] 



as 



/(F 2 ,F 4 ) = F 2 + 2F 2 1/2 F 4 1/3 |A 2 



-1/2 1 > .1/3 



1/2 



After these preparatory moves we can recast our methodology in a convenient specialized 
fashion, suitable for the task at hand. We just face the simple two-equations system: 

I V/(F 2 ,F 4 ) = /iV0(F 2 ,F 4 
{ 0(F 2 ,F 4 ) = 1 



(30) 



where V = (<9f 2 , 9f 4 )- Straightforward solution of it yields 



F 2 " 1/2 (l - F,)- 1 / 3 (7F 2 - 3) = 3 |A 2 | 1/2 |A 4 p 1/3 , F 4 = (1 - F 2 f , (31) 

from which we obtain F 2 and _F 4 . Substituting them into ( 129]) we determine a and, of 
course, the eigenvalue E = a/8. Consider now our SE (128]) . taking k = 1 and a given value 
of A (0.0001 < A < 10000). The function [Cf. Eq. [26] f{F 2 ,F 4 ) = I (x 2 ) exhibits, as a 
function of its arguments, a unique "critical" point that satisfies (131]) . Using / = f critical, 
that optimizes the CR-bound, we find a ground-state eigenvalue that is in good agreement 
the literature. In this way, after properly dealing with (12"^]) . with the Fj. regarded as "FIM 
statistical weights" that optimize the Cramer Rao inequalities, we determine a as a function 
of the Afc without passing first through a Schrodinger equation, which is a notable aspect of 
the present approach. 



Interestingly enough, the Cramer-Rao inequality us equivalent to the quantum uncertainty 
principle (see the Appendix for details and references). Thus, our methodology actually 
employs Heisenberg's celebrated principle to pick up just one solution among the several 
ones that our partial differential equation possesses. 



Table: Ground-state eigenvalues of the SE ()28[) for k = 1 
and several values of the anharmonicity constant A. The 
values of the second column correspond to those one finds 
in the literature, obtained via a numerical approach to the 
SE. These results, in turn, are nicely reproduced by some 
interesting theoretical approaches that, however, need to 
introduce and adjust some empirical constants [121 ]. Our 
values, in the third column, are obtained by means the 
present theoretical, parameter-free procedure. The fourth 
column displays the associated Cramer-Rao bound. 



A 




E = a/8 


f = I{x 2 ) 


0.0001 


1.000074 


1.000074 


1.000059 


0.001 


1.000748 


1.000739 


1.000591 


0.01 


1.007373 


1.007263 


1.005824 


0.1 


1.065285 


1.063047 


1.051255 


1 


1.392351 


1.353533 


1.296590 


10 


2.449174 


2.213973 


2.040974 


100 


4.999417 


4.212932 


3.782394 


1000 


10.639788 


8.587748 


7.599439 



5. CONCLUSIONS 

On the basis of a variational principle based on Fisher's information measure, free of ad- 
justable parameters, we have obtained the Schrodinger energy-eigenvalues for the funda- 
mental state of the quartic anharmonic oscillator (for several anharmoniticy- values). Our 
theoretical results, obtained without passing first through a Schrodinger equation, are in a 
good agreement with those of the literature. This constitutes an illustration of the power of 
information-related tools in analyzing physical problems. 

Thus, we have in this communication introduced a new general technique for eigenvalue- 
problems of linear operators, whose use seems to constitute a promising venue, given the 
results here displayed. 
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FQM-207 of the Junta de Andalucia (Spain, EU). 
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APPENDIX: CRAMER-RAO AND UNCERTAINTY PRINCIPLE 



It is well known that the Cramer- Rao inequality may be regarded as an expression of Heisen- 
berg's Uncertainty Principle (See, for instance, Remember that a precise statement of 
the position-momentum uncertainty principe reads [l6] 

{Ax){Ap)>^ or (Ax) 2 (Ap) 2 >^ (32) 

where 

(Ax) 2 = ((a; - (x)f) = (x 2 ) - (xf (33) 

{Apf = ((p_( J9 )) 2 ) = (^ ) _ (p) 2. (34) 

In a one- dimensional configuration-space, if ijj is a normalizable real wave function, 

d c d ft r d 

^ = (~ i% d^ = ~ i% I ^dx~^ dX = ~ l 2 J Ihf dX = ° ' (35) 

W = (^&} = ^ < 36 > 

Substituting fl35|) and fl36l) in fl34|) and using ([7]) leads to the above mentioned connection 
between the uncertainty in momentum Ap and the Fisher's measure /, i.e., 

(Ap) 2 = -h 2 J^^dx = ^I. (37) 

If this relation is substituted into ( l32l) we immediately arrive to the the CR-bound, 

h 2 

{Ax) 2 (Ap) 2 > -j — > /(Ax) 2 > 1. (38) 

Coming now back to the {a, Afc}-scenario, one easily ascertains that Eq. (I26|) can be given 
a clear "Heisenberg's aspect" 

M 7, 

1 [ p k |A fc | f r2 - L l ( >■ ■ 2 i\ J- ' 2 ; :! i \ .i- 2 

k=l } 



f(F h ■ ■ ■ , F M ) = £ ~ [ F k |A fc | ] 2 ™ (F.^IA.r 1 / 2 - F^IA, 
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